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ABSTRACT 

We present a numerical investigation of the tidal disruption of white dwarfs by 
moderately massive black holes, with particular reference to the centers of dwarf 
galaxies and globular clusters. Special attention is given to the fate of white 
dwarfs of all masses that approach the black hole close enough to be disrupted 
and severely compressed to such extent that explosive nuclear burning can be 
triggered. Consistent modeling of the gas dynamics together with the nuclear 
reactions allows for a realistic determination of the explosive energy release. In 
the most favorable cases, the nuclear energy release may be comparable to that 
of typical type la supernovae. Although the explosion will increase the mass 
fraction escaping on hyperbolic orbits, a good fraction of the debris remains to 
be swallowed by the hole, causing a bright soft X-ray flare lasting for about a 
year. Such transient signatures, if detected, would be a compelling testimony for 
the presence of a moderately mass black hole (below IO^Mq). 



1. Introduction 



White dwarfs, the end point of stellar evolution for stars with masses from about 0.07 
to 10 solar masses, are extremely common (jWeidemannI Il990l ) . And we don't have to look 
far away, eit her. There are several billions of them in the halo of our very own Milky 
Way galaxy (IReidl l2005l ). They are not only observed in isolation but in binary systems, 
with normal stellar companions, and, less frequently, with compact stellar companions. At 



a distance of tens of kiloparsecs, a number of globular clusters ( iBrodie &: Straderl 12006 



Moehler fc Bondl2008l ) have a high enough central density to let white dwar fs interact and 



collide with other stars, or, if deyelope d, with a central massive black hole (iGerssen et al. 



2OO2I . I2OO3I : iGebhardt et al.ll2002l . |2005| ). And at a distance of less than ten kiloparsecs, our 



- 2 - 



galactic nucleus, a central massive black hole of a few million solar masses (iGhez et al.lll998l ). 
is surrounded by swarms of all kind of stars, some of them white dwarfs and all of them prone 
to collisions. In such dense environments, stars of all varieties can exchange mass, disrupt 
each other or merge, and their merger products can get involved in similar interactions. 
Because of this reason, the inner most ~ 0.1 pc of our Galaxy can be considered an efficient 



stellar colhder (lAlexandeii l2005l ). Stellar dynamics is perfectly adequate in modeling the 
motions of stars as point masses moving under the influence of gravity, even in dense stellar 
systems, unless individual stars closely approach each other or get so close to the central 
black hole that they become vulnerable to tidal distortions. Once this happens, the internal 
structure of the stars becomes important for the further dynamical evolution. This paper 
explores the observational manifestations of such phenomena, with particular reference to 
white dwarfs. 

A white dwarf interacting with a neutron star or a black hole cannot be treated as a 
point mass if it gets so close to the compact remnant that it becomes tidally deformed. The 
tidal interaction between a star and a black hole is mainly characterized by three length 
scales: the stellar radius, -Rwd, the gravitational radius of the black hole, Rg = 2G'Mbh/c^ — 
3 X 10^^Mbh,6 cm, where Mbh,6 denotes the hole's mass in units of lO^M©, and the tidal 
radius, Rr- The tidal radius, defined as the distance within which a stars gets disrupted, 
obviously depends on the type of star being considered. For a white dwarf it is roughly given 
by 

lll /\/fl/3 f Rwd \ ( ^ 



~ 1.2 X 10^Wb'h% cm. (1) 

Essentially, Rr is the distance from the hole at which M-qyi/R^ equals the mean density of the 
passing star. is, in order of magnitude, the same as the Roche radius, which is a precisely 
defined quantity, but only strictly applicable to a star in a circular orbit with synchronized 
spin. The strength of the tidal encounter is measured by the dimensionless parameter 

fM^RlV'' 

which is simply the square root of the ratio of the surface gravity of the star to the tidal 
gravity at the surface when the star is at pericenter distance Rp. Most literature on tidal 
disruption uses (3 = Rr/Rp = r/"^/^ to measure the encounter strength. When t] < 1, the 
star is disrupted in a single fiyby. The energy required to tear the star apart (that is, the 
star's self-binding energy) is supplied at the expense of the orbital energy, which at Rr is 
larger by {MBu/M^df/^ 



For white dwarfs, the ratio between the tidal disruption radius and the gravitational 
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Note that this is indeed inside the gravitational radius for black hole masses exceeding 
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(3) 



(4) 



For this reason, white dwarfs only expe rience disruptive physical conditions when approach- 
ing a moderately massive black hole (ILuminet &: PichonI Il989l ). If Mbh "C MBH.Um it is 
however, sufficiently far outside R^ that the disruption can be approximated as a Newtonian 
process, and it makes little difference whether the hole is described by a Schwarzschild or a 
Kerr metric. When the central hole has a mass Mbh ^ M^d, the size of the star remains 
smaller than R^-. A white dwarf cannot thus be disrupted (Figured]) without entering the 
strong relativistic domain if 



j3 > min 
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The type of the black hole (Schwarzchild or Kerr?) then has an important quantitative 
effect, as does (for a Kerr hole) the orientation of the stellar orbit relative to the spin axis. 

Much of our effort in this paper is therefore dedicated to investigate the disruption of 
white dwarfs by mo derately massive bl ack holes, with particul ar reference to the centers of 



dwarf galaxies (e.g. iMathur et al.l 120081 ) . globular clusters (e.g. iPortegies Zwart et al.l 12004 



Baumgardt et al.l 120041) ari d the intermediate-mass black hole candidates in active galactic 
nuclei ( iGreene fc Hoi I2004J ) . Special attention is given to the fate of white dwarfs of all 
masses that approach a black hole close enough to be disrupted and severely compressed 
to such extent that explosive nuclear burning can be triggered. Tidal disruption of white 
dwarfs has previously been difficult to model due to its three-dimensional nature and the del- 
icate interplay between gravity, gas dynamics and nuclear reactions. Cr ucial aspects of the 
physic s of white dwarf disruption and ignition were first underst ood bv iLuminet fc Pichon 
( 1989), mainly bas e d on the affine model earlier introduced by Carter fc Luminetl (jl983l ). 
Luminet fc PichonI (119891 ) suggested that for very close encounters {rj < 1), relativistic ef- 
fects and compression into the orbital plane could trigger significant nuclear energy release 
in the core (some authors, particularly Wilson & Mathews 2004 and Dearborn et al. 2005, 
believe such considerations could also apply to weakly disrupted encounters for which gen- 
eral relativistic effects can overwhelm the stabilizing tidal forces and still cause relativistic 
compression). The affine model of Carter and Luminet approximates the white dwarf as 
a time-dependent ellipsoid, with a fixed density profile. Such a model is of limited utility 
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if the true matter distribution deviates from the assumed shape and density profile. To 
study the tidal disruption problem, we use a three dimensional smoothed particle hydro- 
dynamics method (SPH) to solve the equations of hydrodynamics. This method makes no 
restrictive assumptions on the hydrodynamic degrees of freedom and allows us to study the 
three-dimensional encounter with full generality. Due to its Lagrangian nature SPH is per- 
fectly suited to follow tidal disruption processe s during which the co rresponding geometry, 
densities and time scales are changing violently ( Rosswog et al. 2008a| jbl). What is more, our 



consistent modeling of gas dynamics together with nuclear energy generation allows for a 
realistic determination of the explosive energy release. 

In this paper, a comprehensive approach is adopted: we follow the tidal disruption and 
compression leading up to the ignition of the white dwarf (when triggered); the accompany- 
ing nucleosynthesis and expansion of the debris; and the subsequent accretion of the bound 
debris. Some of the questions at the forefront of attention include the effects of initial com- 
position and mass of the white dwarf, the passage distance, and the black hole's mass on 
ignition and disruption, and the degree to which the post-disruption dynamics is modified 
by nuclear energy release. We address all of these issues here. 

A detailed description of the numerical methods and the initial models is given in §2J De- 
tailed hydrodynamic simulations of the disruption of white dwarfs of various masses, initial 
composition and passage distances are presented in §^ |l]and[5l the role of the mass of the 
black hole in shaping the evolution and ignition of the disrupted white dwarf is discussed in 
^gj The resulting gravitational wave signals are shown in §21 Discussion and conclusions are 
presented in §71 



Numerical Methods and Initial Model 



Th e observational consequences of stellar disruption depend o n what happens to the 
debris ( iReesI Il988l : lEvans fc Kochanekl Il989l : iKobayashi et al.ll2004l ). To quantify this, we 
have performed detailed three-dimensional, hydrodynamical calculations. The gas dynamics 
is coupled with a nuclear network to explore the effects of n uclear energy generatio n dur- 
ing the strong compression phase. The reader is referred to iRosswog et al.l (j2008al ) for a 
complementary description of the numerical methods employed in calculating the disruptive 
event. 
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2.1. Hydrodynamics 



The SPH-formulation used in this study follows c 
a brief derivation of the equations can be foun d in 



■osely the one desc r ibed i n iBenzl ( Il990l ) , 
Rosswog fc Pried (120081 ) . Forces from 



self-gravity are calculated using a binary tree (IBenz et al.l Il990bl ) while the gravitational 
forces from the central bl ack hole are calculated using a Paczyhski-Wi ita pseudo potential 
( jPaczyhski fc Wiitalll980l ) with a polynomial extension (IRosswog 1120051 ) to avoid the singu- 
larity at the Schwarzschild radius. We have taken particula r care to avoid artifacts from the 
use of artificial viscosity (AV). The so-called Balsara-switch (jBalsaralll995l ) is implemented to 
avoid spurious shear forces . More importantly, we use time dependent viscosity parameters 
( iMorris fc MonaghanI 119971 ). so that far from shocks artificial viscosity is essentially absent, 
but near a shock front the associated v iscous parameters rise to values that are able to avoid 
post-shock oscillations (see Fig. 1 of (IRosswog et al.ll2008al ) for an illustration). Once the 
shocks has passe d, the parameters de cay again. Details of the AV implementation and tests 
can be found in (jRosswog et al. 2000). 



A MacCormack predictor-corrector method is used to evolve the fluid in time. The 
method is implemented with individual time steps, i.e. each particle i is advanced on its own 
time step dU = 2"'' x dtmm, where citmin is the smallest step of all the particles, rii is chosen 
to be the largest integer satisfying the condition dti < dti^des, where dti^des is the desired time 
step for particle i. With this time marching implementation the total energy is conserved 
to better than 4 x 10~^ and the total angular momentum to better than 2 x 10^^. Note 
that this could, in principle, be improved even further by tak ing into account the so-called 
"grad-h" -terms (jSpringel fc HernquistI 120021: iMonaehanl l2002l) and extra-terms arising from 
adapting gravitational smoothing terms (jPrice fc MonaghanI 120071 ) . A co mparison between 
the fo rmulation used here and the one using "grad-h" -terms can be found in IRosswog &: Price 
(|2007h . 



2.2. Equation of State 

We use the HELMHOLTZ equation of state (EOS), developed by the Center for As- 
trophysical Thermonuclear Flashes at the University of Chicago. The EOS allows to freely 
specify the chemical composition of the gas and it can be coupled to a nuclear reaction 
network. The electron-/positron equation of state is calculated without approximation. 
In other words, it makes no assumption about the degree of degeneracy or relativity and 
the exact expressions are integrated numerically. The nuclei in the gas are treated as a 
Maxwell-Boltzmann gas, the photons as blackbody radiation. The EOS is used in tabular 
form with densities in the range 10^^° < pYg < 10^^ g cm~^ and temperatures varying be- 
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tween 10^ and 10^^ K. A sophisticated, biquintic Hermite polynomial interpolation is used to 
enforce the thermodynam ic consistency (i.e. the Maxwell-relations) at interpolated values 
Jximmes fc Swestv lEoOof ). 



2.3. Nuclear Burning 



To account for the feedback onto the fluid from nuclear transmutations we use a minimal 



nuclear reaction network developed by lHix et al.l (119981 ). It couples a conventional a-network 
stretching from He to Si with a quasi-equilibrium-reduced a-network. The QSE-reduced 
network neglects reactions within small equilibrium groups that form at temperatures above 
3.5 X 10^ K to reduce the number of abundance variables needed. Although a set of only 
seven nuclear species is used, this network reprod uces all burning stages from He-burning to 
NSE accurately. For details and tests we refer to iHix et al.l (119981 ). 



The network is coupled to the hydrodynamics in an operator splitting fashion, i.e. 
hydrodynamics and nuclear burning -which may in extreme cases require vastly different time 
steps- are integrated separately. In a first step, the hydrodynamic equations are integrated 
with the above described predictor-corrector scheme to obtain new quantities at time step 
fn+i^ In this step we ignore the nuclear source term in the energy equation, the new value for 
the specific energy of particle a is denoted by u^"*"^. The specific energy has to be corrected 
for the energy release that occurred between and t""*"^: 



-a,n— >nH-l 



m + l 



dt 



{pa{t),Ta{t),Yk,ait)) dt 



(6) 
(7) 



where A'a is Avogadro's constant, rrijC^ is the mass energy of nucleus j and Yj^a is the 
abundance of nucleus j in particle a. We use Pa(t) ~ Pa(t"') + {Paif^^^) — pa{t^)} 



and Ta{t) 



t-t" 



-{Ta{t 



n+l^ 



T(j(t")}, where Ta{t^~^^) is the temperature derived from 



m""*"^, to integrate the abundances K, „ via the impli cit backward Euler method (the network 



integration is described in detail in (iHix et al. The final value for the specific energy 

at time t""*"^ is given by 



u 



n+l 



~n+l _|_ , 



^Note that the tempe rature evolution along a hydrodynamical time step is different from the description 
in iRosswQg et al.l (|2008al) . In practice, we only see tiny differences between both approaches. 
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The EOS is then called again to make all thermodynamic quantities consistent with this new 
value m""*"^. Once the derivatives have been updated, the procedure can be r epeated for the 



next ti me step. For details of the time step criteria we refer the reader to iRosswog et al. 
fl2008ah . 

As will be shown below, the peak compression occurs at a spatially fixed point (see density 
panels in Fig. [6]). The white dwarf fuel is fed with free-fall velocity vs = {2GMbk/ RpY^'^ = 
1.6 ■ 10^ km s~^ (Mbh/1000 Mq)^''^ (i?p/10^ km) ^^'^ into this compression point. This is 
many orders of magnitude larger than typical fiame propagation speeds, therefore fiame 
propagation effects can be safely neglected for this investigation. 



2.4. Initial Model 

It is a vital ingredient for any calculation to start out from initial conditions that are 
as accurate as possible. In order to build stars in hydrostatic equilibrium, it is therefore 
necessary to find SPH particle distributions whose number density refiects the equilibrium 
mass density distribution. To this end we solve the spherically symmetric Lane-Emden 
equations for a star of given mass and composition to find a one dimensional density profile. 
For simplicity, and to start from a conservatively low value, we assume a uniform temperature 
of To = 5 X 10"^ K. This assumption will be relaxed in one run (run 10, see Table 1) to test 
for the sensitivity to this initial condition. Once the solution to the Lane-Emden equation 
has been obtained, we distribute the desired number of particles, N, inside a unit sphere 
according to a close-packed configuration. Subsequently, we map this distribution into the 
volume of the star that is to be constructed. As an example, the particle distributions before 
and after mapping of a 0.2 Mq white dwarf are shown in Figure El 

The mapping of the initial configuration is done so that the SPH-particle number density 
in the star, n, yields 

pLE(r) = m ■ n{r), (9) 

where m is the mass of each particle and ple is the density profile obtained by solving the 
Lane-Emden equation. We denote quantities referring to the star (unit sphere) with (un- 
) primed variables. If the unit sphere has a constant density po and contains particles, each 
with mass m = M^^/N, the mass in the shell between r„_i and is M„ — 47rr^_]^Arpo- As 
illustrated in Figure [31 the image of this shell (in the star) contains the same mass, but now 
between r^_^ and r^, so that ^ 4:nr'^_-^Ar'^pn-i. The condition = M„ then yields 



(10) 
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The resultant configuration is very close to hydrostatic equihbrium. To find the "real" 
numerical equilibrium state we relax this configuration with the full hydrody namics code by 



applyi ng an artificial, velocity dependent damping force as, for example, in iRosswog et al. 



too4 \ 



We use this procedure for all but the heaviest white dwarfs (1.2 Mq). Since their 
adiabatic exponent is already approaching the critical 4/3, these stars are very centrally 
condensed and, as a result, the bulk of the SPH particles resides very close to each other in 
the stellar center. This would, even for a moderate resolution of a few hundred thousand 
particles, result in most of the particles running at time steps that are orders of magnitude 
smaller than the dynamical time of the star. To alleviate this computational burden in the 
heaviest white dwarfs we allow for slightly different particle masses (about a factor of 10 
from the center to the stellar surface). 

Throughout this simulation set we assume a uniform nuclear composition accross the white 
dwarfs. Stars with masses < 0.6 Mq are instantiated as pure helium, more massive stars 
are modeled as 50% carbon and 50% oxygen. Note that helium core burning may produce 
composition gradients throughout the stars with oxygen being more abundant in the stellar 
core. Since it is mainly the core that becomes ignited during the compression, see Sect. El 
such a gradient may slightly reduce the nuclear energy that can be released in principle. 
Common to all calculations is the initiation of the calculations with the white dwarf being 
placed safely outside Rr. Initial separations are at least twice, in most cases three times 
the tidal radius to allow the initially spherical star to adjust properly to the changing tidal 
potential as it approaches the black hole. All stars are set onto parabolic orbits so that Rp 
is larger than Rg. 

The performed runs are summarized in Table [H 



3. Events in the Life of a Tidally Disrupted White Dwarf 

Since the focus of this study is closely related to the nuclear energy release during strong 
encounters, practically all following calculations refer to cases with penetration factors f3 
clearly beyond unity, see Table 1. To illustrate the evolution of a white dwarf that is just 
marginally disrupted, we show in Fig. H] a 0.6 Mq CO white dwarf that passes a 1000 Mq 
black hole with a penetration factor of f3 = 0.9 (run 10, Table 1; keep in mind that the tidal 
radius as defined above is approximate and its exact value depends in the internal structure 
of the disrupted star). The upper panel shows snapshots at t=0.34, 3.43, 6.86, 10.29, 13.72, 
17.15, 20.58 and 24.01 s after the simulation start at a separation of 3 Rr- The lower panel 
shows late stages at t= 138.24 and 380.53 s. The circle indicates in both cases the location 
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of the tidal radius. At about 1.5 Rr the star becomes noticeably deformed and substantially 
spun up by the time it reaches pericentre (snapshot 5, upper panel). While receeding from 
the black hole large, puffed-up lobes form at the extremes of the star. The inner lobe returns 
to the black hole (see t=138.24 s, lower panel) while the outer is ejected to infinity, both are 
connected by a well-defined and homogenous spaghetti-like tube of white dwarf debris. 
The opposite limit, an exremely strong encounter, is illustrated in Figure [51 The snapshots 
are taken from our numerical simulations of a 0.2 Mq white dwarf approaching a IO^Mq black 
hole on a parabolic orbit with pericenter distance well within the tidal radius (i?p = Rr /12; 
run 1 in Table 1). While falling inwards towards the hole, the white dwarf develops a 
quadrupole distortion, which attains an amplitude of order unity by the time of disruption 
at -R ~ Rr- 

When the star is still far from the black hole, it remains close to its initial stationary 
equilibrium state as characterized by the usual virial relation. This is because the timescale 
characterizing the rate of change of the tidal force will initially be very long compared with 
the intrinsic timescale characterizing the corresponding quadrupole oscillations of the star. 
The tidal bulge raised on the star by the black hole becomes an order unity distortion 
near the tidal radius. The resultant gravitational torque spins it up to a good fraction of 
its corotation angular velocity by the time it gets disrupted. The large surface velocities 
and the order unity tidal bulge combine to overcome the star's self gravity and lead to the 
disruption of the star. 

Our present investigation is concerned with cases of deeper penetration within the 
tidal radius so that the core as well as the envelope are affected. The behavior of a white 



dwarf passing well within the tidal radius exhibits special features (ICarter fc Luminetlll983 



Bicknell fc Gingoldl 119831 ). As illustrated in Figure [6], the degenerate star is not only elon- 
gated along the orbital direction but also severely compressed perpendicular to the orbital 
plane. This anisotropy can be understood as arising from the fact that the principal tidal 
axes within the orbital plane will rotate prior to pericenter passage so that the corresponding 
effects of elongation and compression will roughly cancel out, whereas the third principal axis 
retains a fixed direction so that compression orthogo nal to the orbital plane is uncontested 



see, e.g., the appendix of (IBrassart fc Luminetll2008l )). Each section of the star is squeezed 



through a point of maximum compression at a fixed point on the star's orbit. This takes 
place on a timescale comparable to the crossing time of the star through periastron, 

r. Rwd no/ \ f Rwd f ^BH \ /-.-.N 



V, 
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where 

^ iW-c . 5 X 10» (^) (^) (^) ,12) 

is the orbital velocity at periastron. 

During this very short lived phase, the star attains its maximum degree of compression 
(here by a factor ~ 100) although its orbital-plane section increases only by a factor of a few 
due to the continually changing directions of the two principal tidal axis (Figure E])- The 
distortions A_R„d/-Rwd ~ 1 impart supersonic bulk flow velocities of order 



during t he drastic compression of the stellar material. Th is compression is thus halted by 



a shock (IKobayashi et al.ll2004l : iBrassart fc Luminetl 120081 ). raising the matter, which then 



rebounds perpendicular to the orbital plane, to a higher adiabat (Figure [7]). As a result, 
the temperatures increase sharply and trigger explosive burning (of He for the case shown 
in Figure [6]). 

The typical temperature of the shocked star and the thermal energy produced by shock 
heating can be roughly estimated from the virial theorem as 

r~6xlO»f-^V7j^U, (14) 



lO^cm ) \ 0.6 Mq 



and 



2 



^..™-^c^lO-(^)" (jf^j e.g. (15) 

Adiabatic compression alone can only increase the stellar surface and interior layer temper- 
atures by a modest factor. If adiabatic compression was the only source of heating, the 
response to the flow to the varying potential would be ~ Rwd/cg ~ 4 s, which is significantly 
longer than the dynamic crossing time. 

The evolution of the abundances (for the simulation shown in Figure [5l) is shown in Fig- 
ure [HI within the 12 s during which the entropy increases the fastest. During the about one 
second long period of compression, the temperature increases beyond 3 x lO^K, approach- 
ing but not quite reaching nuclear statistical equilibrium (NSE). During this brief period 
of compression, nuclei up to and beyond Si are synthesized. The initial composition was 
pure He and the final mass fraction in iron-group nuclei is about 15%. This result should be 
taken as a modest underestimate, since the seven species nuclear network only provides an 
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approximation to the detailed nuclear processes. Post-processing calculations, using a 300 
isotope nuclear network over thermodynamic particle histories resulting from these calcula- 
tions, show significant nuclear fiow above silicon, for helium-rich portions of the gas with 
peak temperatures above 2 x lO^K. As a result, heavier elements (like calcium, titanium and 
chromium) would be made, accompanied by a modest increase in the energy generation. It 
is thus safe to conclude that the white dwarf is tidally ignited and that a sizable mass of 
iron-group nuclei is injected into the outfiow. 

The variation of the specific energy in the released gas, in the absence of explosive energy 
input, is determined mainly by the relati ve depth o f a mass element across the disrupted 



star in the potential well of the black hole (jReeslll988l ). This is much larger than the binding 



energy and the kinetic energy generated by spin-up near pericenter. Even though the mean 
specific binding energy of the debris to the hole is comparable with the self-binding energy 
of the original star, the spread about this mean is larger by a factor {M-qyi/ M^^y/'^ . Nuclear 
energy released during the drastic compression and distortion of the stellar material, further 
enhances the spread in specific energies at pericenter (Figure [9]), and, as a result, the mass 
fraction escaping on hyperbolic orbits is increased from ~ 50% to ~ 65% of the initial mass 
of the white dwarf. The mass fraction that is ejected rather than swallowe d, though less 



spectacular than typical Type la supernovae (iHillebrandt fc Niemeyerl |2000| ). should have 
many distinctive observational signatures (the reader is referred to Kasen et al. 2008 for 
a detailed description of the optical light curves and spectra resulting from the unbound 
debris before it becomes translucent). First, the explosion itself should be different, since 
the disrupted, degenerate stars should be, on average, lighter than those exploding as type 
la supernovae. Second, the spectra should exhibit large Doppler shifts, as the ejected debris 
would be expelled with speeds > 10^ km/s. Finally , the optical light curve should be rather 
unique as a result of the radiating material being highly squeezed into the orbital plane (one 
thus expects different timescales for conversion of nuclear energy to observable luminosity 
when compared with normal type la events). 

Although the explosion will increase the fraction of ejected debris, a good fraction (~ 
35%) remains to be accreted on to the hole. The returning gas does not immediately produce 
a flare of activity from the black hole. First material must enter quasi circular orbits and 
form an accretion torus faees 1988 : Evans fc Kochanek 1989). Only then will viscous effects 



release enough binding energy to power a flare. The bound orbits are very eccentric, and the 
range of orbital periods is large. For white dwarfs, the orbital semi-major axis of the most 
tightly bound debris is 

^BH A ^^'^ ( Rwd \ f ^ '^^^ 
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and the period is only 

\300RJ \10^MqJ ^ ' 

If the gaseous debris suffered no internal dissipation due to high viscosity or shocks, it 
would, after one or two orbital periods, form a highly elliptical disc with a big spread in 
apocenter distances between the most and least bound orbits, but where at peri center, Rp, 
the radial focusing of the orbits acts as an effective nozzle (Figure HI and [TOj) . After pericenter 
passage, the outflowing gas is on orbits which collide with the infalling stream near the 
original orbital plane at apocenter, giving rise to an angular momentum redistributing shock 
(Figure [TOl) much like those in cataclysmic variable systems. The debris raining down would, 
after little more than its free-fall time, settle into a disc. This orbiting debris starts to form 
when the most tightly bound debris falls back. The simulation shows that the first material 
returns at a time < ta? with an infall rate of about IO^Mq yr^^ (Figure [TT]) . Such high 
infall rates are expected to persist, relative steadily, for at least a few orbital periods, before 
all the highly bound material rains down. The vicinity of the hole would thereafter be fed 
solely by i njection of the infall i ng matter. The early mass infal l rate is sensitive to the stellar 



structure (ILodato et al. 



2008 



Ramirez- Ruiz fc Rosswod l2008l ) . at late times, t > tfb ~ 600 



s, it drops off as t-^/^ jR,eesl l llQSsI : |Phinneylll989l ). Once the torus is formed, it will evolve 



under the influence of viscosity, radiatively cooling winds and time dependent mass inflow. 

A luminosity comparable to the Eddington value, ~ i^Edd = 10^^(Mbh/10^Mq) erg s~^, 
can therefore only be maintained for at most a year; thereafter the flare would rapidly fade. 
It is clear that most of the debris would be fed to the hole far more rapidly than it could 
be accepted if the radiative efficiency were high; much of the bound debris must either 
escape in a radiatively-driven outflow or be swallowed inefficiently. The rise and the peak 
bolometric luminosity can be predicted with some confldence. However, the effective surface 
temperature (and thus the fraction of luminosity that emerges predominantly in the soft 
X-ray band) is harder to predict, as it depends on the size of the effective photosphere that 
shrouds the hole. Such transient signals, if detected, would be a compelling testimony for 
the presence of moderately massive black holes in the centers of globular clusters and dwarf 
galaxies. 



4. Influence of Black Hole's Mass on Disruption 



The characteristic tidal radius, Rr, for a given white dwarf star is solely determined 
by the black hole's mass, while the strength of the tidal encounter is traditionally measured 
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by the dimensionless parameter (3 = Rr/Rp (provided (3 < (3g). The aim of this section is 
to illustrate, with the help of a few specific calculations, the role of the black hole's mass 
in shaping the evolution and ignition of the disrupted white dwarf. To quantify this, we 
have performed detailed three-dimensional, hydro dynamical calculations of the dynamics of 
a 0.2Mq white dwarf approaching black holes of various masses on parabolic orbits with 
p = 5. For black hole masses IO^Mq < Mbh < IO^Mq and p = 5, a 0.2Mq white dwarf 
can be disrupted without entering the strong relativistic domain (Figure [Tj). However, com- 
plete disruption occurs sufficiently close to the hole for a Newtonian approximation to be 
inadequate (Figure [T2|) . 

The behavior of a 0.2Mq white dwarf approaching black holes of various masses is 

illustrated in Figure [T31 As discussed previously, the effects of the black hole's mass (for 

1 /3 

a fixed P and M„d) will be to move the critical pericentric distance by Mgjj and change 

— 1/3 

the star's crossing time through the maximum compression point by Mg^ . As a result, 
stars approaching black holes of increasing mass will appear more elongated along the orbital 
plane, as each section of the star is squeezed faster through a point of maximum compression 
at a fixed point on the star's increasingly extended orbit. Such stars are not only more 
elongated but are even more severely compressed into a prolate shape. If the initial white 
dwarf is made of pure He (as in the 0.2 Mq case), the co mbustion rate will be dete rmined 



by the 3a reaction, on time scale approximately given by (iKhokhlov Sz Ergmalll986l ) 

tb ^ 9.0 X IQ-'T^ exp(4.4/T9) pf, (18) 

where Tg is the temperature in units 10^ K and pe the density in 10^ gcm~^ . 

If the time scale on which the white dwarf can react, its dynamical time scale tQ = 
(Cpwd)"^^^, is much shorter than the burning time scale tb, the star can expand rapidly 
enough to quench burning by a reduction of temperature and density. Appreciable burning 
will therefore only take place if tb ■'^ ^g- As illustrated in Figure [HI this comparison of time 
scales can be used as a simple estimate for whether or not substantial He-burning will occur. 
Promising models for thermonuclear explosions are those that reach high temperatures at 
high densities (upper right corner). For a fixed P, increasing Mbh raises the nuclear energy 
release. For black hole masses Mbh > 10^ M©, nuclear energy released during the drastic 
compression and distortion of a 0.2 Mq white dwarf approaching moderately massive black 
holes with P = 5 is fast enough to release energy in excess of that required to tear the star 
apart: Eq (Figure [T^ . When the black hole mass is ^ lO^M© most white dwarfs would be, 
however, swallowed whole {Rr < Rg). 
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5. Tidal Disruption for a Variety of White Dwarfs 

White dwarf s are thought to b e the final evolutionary state of all stars whose mass is 



below ~ 10 Mq (IWeidemannI 119901 ) . Stars with masses less than about half a solar mass 
become degenerate before helium ignition, and therefore will end their lifes as helium white 
dwarfs. In isolation, such low mass stars have life times much longer than the present age of 
the Universe. Still, such He white dwarfs are observe d, they are thought to be the result of the 



evolution of a close binary system (e.g. iHeberj |2002|) . Stars of low or medium mass, achieve 



helium burning and become carbon-oxygen white dwarfs (perhaps the com monest sort). 



The o bserved white dwarf mass distribution is strongly peaked around 0.6 Mq (IKepler et al. 



2003). 

When a white dwarf is subject to strong tidal compression, triggering of nuclear processes 
in the stellar cores depends on initial composition. To investigate this dependence, we have 
performed calculations with various white dwarf masses/compositions. As outlined above, 
our white dwarf models are initialized with pure He-composition for M^^ < 0.6 Mq and 50% 
carbon and 50% oxygen otherwise. 

Although the thermodynamical evolution and the nuclear energy release are very sensi- 
tive to the initial stellar composition (in this text, we have so far considered only a pure ^He 
0.2 Mq white dwarf), we found that ignition is in fact a natural outcome for white dwarfs of 
all masses passing well with in the tidal radius. For example, a C/0 I.2M0 white dwarf (here 
assumed to be composed of 50% carbon and 50% oxygen throughout the star) approaching a 
500 Mq black hole on a parabolic orbit with pericenter distance (3 = 3.2 ignites (Figure [T6l) 
and, as a result, at least 0.6 Mq of Fe are synthesised in the fiow (Figure [T7|) . Similarly, at 
least 0.66 Mq of Fe are synthesised when a C/0 I.2M0 white dwarf approaches a IO^Mq 
black hole with /3 = 3. By contrast with the previous two cases, explosive energy release is 
too slow to release much energy on a dynamical timescale when a C/0 I.2M0 white dwarf 
approaches a 10'^ Mq black hole with f3 = 1.5 and, as a result, only 1.4 x 1O~^M0 of Fe are 
synthesised. 

0.18 Mq of Fe are synthesised when a typical 0.6 Mq C/0 white dwarf approaches 
a 500 Mq black hole with (3 = 5. To explore how sensitive our models are to the initial 
temperature, we have re-run this simulation, but now with a hot white dwarf. The initial 
specific energy of a cold (T = 5 x 10^ K) equilibrium model was adjusted in a way that 
the temperature varied linearly between 5 x 10'' K in the centre to 5 x 10^ K at the stellar 
surface. In the subsequent relaxation the star expanded slightly, since the degeneracy was 
lifted to a small extent. As a result of the reduced degeneracy, lower peak temperatures 
and densities are reached (see Figure [TSll. reducing slightly the amount of nucleosynthesis 
which takes place. We consider the shown results as representative, but the actual range of 
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behaviour for given sets of masses and orbits may be wider due to gradients in temperature 
and nuclear composition which have not been explored systematically here. To summarize, 
in the most favorable cases, the nuclear energy release, is comparable to that of typical Type 
la supernovae. 

In the tidal pinching process, explosive nucleosynthesis is likely to proceed only for 
white dwarfs passing well within the tidal radius (Figure [19]). Explosive energy release, 
as calculated here, appears to be a natural outcome for 0.6 Mq (1.2 Mq) white dwarfs 
approaching moderately massive black holes with f3 > 5 (/? > 3). 



6. Gravitational Waves 



A white dwarf approaching the tidal radius will be disrupted in a single flyby. The 
resulting stellar debris trail is not compact enough to emit strong gravitational waves after 
leaving Rp. The detectable gravitational wave signal will therefore have a burst-like behavior, 
roughly characterized by an amplitude h and a duration t ~ 1//, where 



h 



and 



GQ 
c^D 



c^pD 



7 X 10-2° (3 



D 



lOkpc 



-1 



wd 



0.6Mq 



4/3 



R 



wd 



109 cm 



2/3 
BH,3' 



/ 



BH 



1/2 



0.3 /?^/2 



wd 



1/2 



R 



wd 



-3/2 



Hz. 



(19) 



(20) 



O.QMqJ VlO^cm^ 

LISA will be able to detect gravitational waves of amplitu de h ~ 10"^^ for burst sources in 



the frequency range / ~ 10 ^ — 10 ^Hz (IDanzmanru l2003l ) . Gravitational waves from white 
dwarf stellar disruption could thus be detectable if M„d > 0.4 Mq and the source distance 
D < 10 kpc (Figure [20]). 



The gravitational wave amplitudes and h-^ shown in Figure [21] are calculated in the 
quadrupole approximation. The reduced quadrupole mo ments can be written in terms of 
the SPH-particle properties (jCentrella fc McMillan] Il993l ) 



(21) 



The second time derivatives, Ijk, can then be expressed in terms of the particle properties 
by simple, direct differentiation. The retarded gravitational wave amplitudes for a distant 
observer along the z -axis at distance D are given by 



Dh. 



'to; 



(22) 
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and 

D/ix=2^4,. (23) 

Both duration and amplitudes of the gravitational wave bursts that are produced by the 
disruption calculations (Figure are in agreement with the simple estimates given above. 



Discussion 



This paper presents a computational investigation of the mechanical and nuclear evo- 
lution of white dwarfs passing well within the tidal radius of a moderately massive black 
hole. A comprehensive top to bottom approach is adopted: we follow the tidal disruption 
and compression leading up to the ignition of the white dwarf; the complex propagation of 
the nuclear energy release through the star; the resultant gravitational wave signal and the 
subsequent accretion of the bound debris. 

This paper has outlined several potentiall y observable effects. The detection of a p ecu- 
liar, underluminous thermonuclear explosion (jRosswog et al.l l2008bl : iKasen et al.l l2008l ) ac- 
companied by a thermal transient signal of predominantly soft X-rays with a peak luminosity 
L ~ LEdd = 10^^Mb//^3 erg/s, fading within a year would, if detected, be a compelling testi- 
mony for the existence of a new mechanism by which white dwarf ignition can be achieved. 
Although the thermodynamical and nuclear energy of the star is sensitive to the initial com- 
position, we found that thermonuclear ignition is a natural outcome for white dwarfs of all 
masses passing well within the tidal radius, with lighter stars requiring deeper penetration 
into the tidal radius due to their lower densities. For simplicity, we have instantiated our 
initial carbon oxygen white dwarf models (Mwd > 0.6 M©) as homogeneously mixed stars 
with a 50% mass fraction of each nucleus. While such internal cher nical profiles are likely 
accur ately realized in nature in very massive white dwarfs (~ 1 M©) (IMazzitelli fc Dantona 
1986I ). for lower masses the gravothermal adjustment of the interior during the cooling phase 
produces oxygen-enhanced stellar cores surrounded by very carbon-rich mantles [Xc ~ 0.8). 
The exact radial distribution depends on the exact value of ^"^Cja, 'jY^O rat e and the details 



of ho w convection proceeds (IMazzitelli &: Dantonalll986l : ISalaris et al.l 119971 : IStraniero et al. 



20031 ). but this general stratification tendency is well-established. Thus, the disruption of a 
standard 0.6 Mq white dwarf should produce a highly carbon-enriched remnant atmosphere. 
Maybe , the recently detecte d carbon-rich transient SCP 06F6 accompanied by an X-ray 
signal (iGaensicke et al.ll2008l ) is already the first example for this class of object. 



The gravitational forces from the central black hole are currently calculated using a 
Paczyhski-Wiita pseudo potential but our goal is to incorporate the effects of general rel- 
ativity, as white dwarfs passing well within the tidal radius cannot be disrupted without 
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entering the strong relativistic regime. The form of the black hole (Schwarzschild or Kerr?) 
then has an important quantitative effect, as does (for a rotating Kerr hole) the orientation 
of the stellar orbit relative to the bl ack hole spin axis. The o rbits are then not ellipses, but 



may turn through 27r or even more (ILuminet fc Pichonlll989l ). A fluid element, which in the 



case of an elliptical orbit would cross the orbital plane just once, may then have two or more 
traversals. This opens up the possibility of multiple shocks. 

For a white dwarf which does not pass close enough to the hole to release much energy on 
a dynamical timescale, adiabatic cooling would severely reduce the internal radiative content 
before the debris became translucent (just as a supernova would be optically inconspicuous 
in the absence of continuing energy injection in the months after the explosion). There would 
be no transient until, as discussed above, the bound debris fell back onto the hole after tfb. 
The integrated output from this flare could, in principle, amount to a few per cent of the 
white dwarf's rest mass, but would probably be significantly less, because most of the debris 
would be fed to the black hole far more rapidly than it could be accepted if the radiative 
efficiency were high; much would then be swallowed inefficiently or most likely escape in 
radiatively-driven directed outflow: its ram pressure and subsequent heating could inhibit 
the steady accretion that would otherwise be inevitable in any galaxy or globular cluster 
harboring a moderately massive black hole. 
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Table 1: Summary of the performed runs. M„d and Mbh are the masses of the white dwarf 
and the black hole, respectively, (3 is the ratio of tidal radius and pericentre distance. The 
type of gravity is indicated by N (Newtonian) or PW (Paczyhski-Wiita). Column six states 
the number of SPH particles used in the simulation, £^burn is the energy generated in burning 
processes. To put this number into context we briefly state the gravitational binding energies 
(in ergs) of the different white dwarfs: log(£'bin,o.2M0) = 49.13, log(i?bin,o.6M,;,) = 50.12, 
log(£^bin,i.2M0) = 50.98. "Fe" labels the mass in iron-group elements, "expl." in the comment 
column indicates that the produced nuclear energy exceeds the WD gravitational binding 
energy. 



run 




Mbh 




grav 


SPH part. 


log(^burn) 


"Fe" [Mq] 


comments 


1 


0.2 


1000 


12 


N 


4034050 


50.46 


0.025 


expl. 


2 


0.2 


1000 


12 


PW 


4034050 


50.44 


0.034 


expl. 


3 


0.2 


1000 


12 


PW 


200452 


50.44 




r = 5/3-polytrope 


4 


0.2 


100 


5 


PW 


100027 


49.57 


< 10-1° 


explore BH influence, expl. 


5 


0.2 


500 


5 


PW 


100027 


49.64 


< 10-10 


explore BH influence, expl. 


6 


0.2 


1000 


5 


PW 


100027 


49.76 


< 10-10 


explore BH influence, expl. 


7 


0.2 


5000 


5 


PW 


100027 


49.93 


< 10-10 


explore BH influence, expl. 


8 


0.6 


500 


5 


N 


502479 


50.68 


0.18 


expl. 


9 


0.6 


500 


5 


N 


502479 


50.62 


0.13 


hot, initial WD 


10 


0.6 


1000 


0.9 


N 


1006446 


0.00 


0. 


no nuclear burning 


11 


0.6 


1000 


5 


PW 


502479 


50.43 


3 X 10-^ 




12 


0.6 


10000 


1.5 


PW 


502479 


45.07 


< 10-10 




13 


1.2 


100 


3.5 


N 


100027 


51.01 


0.58 


expl. 


14 


1.2 


500 


2.6 


PW 


502479 


51.16 


0.66 


expl. 


15 


1.2 


1000 


1.5 


PW 


502479 


49.63 


0.014 




16 


1.2 


1000 


3.0 


N 


502479 


51.10 


0.63 


expl. 
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Fig. 1. — This diagrams depicts the relevant domain where the black hole's tidal effect can 
be disruptive to a white dwarf star. The penetration factor (5 is plotted as a function of 
the black hole mass Mbh- Note that a black hole of ^ IO^Mq can swallow all white dwarfs 
without first disrupting them (adapted from Luminet & Pichon 1989). 
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Fig. 2. — Building initial conditions by solving the spherically symmetric Lanc-Emden equa- 
tions for a star of given mass and composition to find a one dimensional density profile. The 
particles are then distributed inside a unit sphere according to a close-packed prescription. 
Left Panel: Close-packed particle distributions in a unit sphere. This distribution is then 
mapped into the volume of the star that is to be constructed. Right Panel: Particle dis- 
tribution after mapping onto the density profile of a 0.2 Mq white dwarf. Shown are the 
projections of the particle positions onto the xy-plane. This particle distribution is then 
finlly 'relaxed' into its true numerical equilibrium. 
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Fig. 3. — Mapping between unit sphere and star. We denote quantities referring to the star 
(unit sphere) with (un-)primed variables. 




Fig. 4. — Evolution of a 0.6 white dwarf passing a 1000 Mq black hole with a penetration 
factor f3 of only 0.9. The snapshots show t= 0.34, 3.43, 6.86, 10.29, 13.72, 17.15, 20.58 and 
24.01 s after the simulation start in the upper, and 138.24 and 380.53 s in the lower panel. 
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Fig. 5. — A 0.2 Mq white dwarf (modeled with more than 4x 10^ SPH particles) approaching 
a 10^ Mq black hole on a parabolic orbit with pericenter distance Rp = Rr/12 is distorted, 
spun up during infall and then tidally disrupted. Shown are density cuts at various instants 
along the orbital (xy-) plane. Color bar gives the amplitude of logp = [logpmax, logpmin] in 
cgs units: [4.0, 5.3] for t=6.2 s, [3.6, 5.4] for t=6.8 s, and [0, 3.4] for t=10.3 s. 
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Fig. 6. — Tidal deformation of the white dwarf before and after passage through pericenter, 
as the star attains its maximum degree of compression (same simulation as referred to in 
Figure [5]). The panels in the upper rows show cuts of temperature (in units of 10^ K) and 
density (in cgs units) through the orbital (xy-) plane. Color bar gives the values of T = 
[T„,in,T^ax] (logp = [logp^in,logPn,ax]): [0,330] ([3.6,5.4]) for t=Q.8 s, [0,1280] ([1.0,6.0]) 
for t=7.0 s, and [0,3200] ([1.0,6.0]) for t=7.2 s. The panels in the lower two rows show 
the temperature (in units of 10^ K) and column density (in cgs units) distributions in the 
xz-pane (averaged along the y-direction). Color bar gives the values of T = [Tmin,Tmax] 
(logs = [logS^i„,logS^ax]): [0,220] ([12.0,15.7]) for t=6.8 s, [0,720] ([12.0,16.0]) for t=7.0 
s, and [0, 1700] ([12.0, 16.5]) for t=7.2 s. The dimension of the bar scale is 10^ cm. 
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Fig. 7. — Evolution of the entropy e of the central portion of the disrupted white dwarf 
just after it attains its maximum degree of compression (same simulation as referred to in 
Figures [5] and [6]). This compression is halted by a shock, raising the matter to a higher 
adiabat. The panels show entropy (in cgs units) cuts through the orbital (xy-) plane. 
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Fig. 8. — Evolution of the abundances during the disruption and ignition of a 0.2 Mq white 
dwarf passing a IO^Mq black hole (same simulation as referred to in Figures [5H7]) . 
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Fig. 9. — Differential mass distributions in specific energy for tlie 0.2 Mq wliite dwarf debris 
(same simulation as referred to in Figures [5H8]) . 
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Fig. 10. — Column density (in cgs units) and temperature (in units of 10^ K) distributions 
in the orbital plane of the bound debris minutes after disruption. The most tightly bound 
debris would transverse an elliptical orbit with major axis ~ 300i?g before returning to 
Rk, (same simulation as referred to in Figures [5H9]) . 
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Fig. 11. — The rate at which the 0.2 Mq white dwarf debris returns to the vicinity of the 
black hole (same simulation as referred to in Figures IS HTOl) . 
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Fig. 12. — The importance of relativistic effects from the central black hole in determining 
the behavior of white dwarfs passing within the tidal radius. The trajectory of the center of 
mass of a 0.2 Mq white dwarf approaching a 10^ Mq black hole on a parabolic orbit with 
impact parameter (3 are shown using Newtonian gravity (N; black) and the pseudo-Newtonian 
relativistic potential of Paczyhski-Wiita (PW, blue). Each time, the star is modeled with 
more than 10^ partices. 
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Fig. 13. — A 0.2 Mq white dwarf approaching black holes of various masses on a parabolic 
orbit with /3 = 5. The panels in the show density (in cgs units) cuts through the orbital 
(xy-) plane. Color bar gives the values of logp = [logpmin, logpmax]- For Mbh = lO^Mo, 
logp = [4.25,5.3] for t= 14.8 s, [4.0,5.4] for t= 15.4, and [3.75,5.8] for t= 16.1. For Mbh = 
5 X lO^M©, logp = [4.15,5.3] for t= 15.6 s, [3.9,5.65] for t= 16.4, and [3.75,5.8] for t= 17.3. 
For Mbh = lO^M©, logp = [4.1,5.4] for t= 15.6 s, [3.9,5.65] for t= 16.3, and [3.4,5.6] for 
t= 17.2. For Mbh = 5 x lO^M©, logp = [4.0,5.4] for t= 15.4 s, [3.7,5.75] for t= 16.0, and 
[3.4,5.8] for t= 16.3. 
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Fig. 14. — The importance of the central black hole's mass in determining the behavior of 
a 0.2 Mq white dwarf passing within /3 = 5 (same simulations as referred to in Figure [T3|) . 
The evolution of the compressed, and tidally disrupted white dwarf in the p — T plane. The 
hottest 10% of the particles are identified and their average temperature (in units of 10^ K) 
is plotted as a function of their average density (in cgs units). These trajectories always 
start cold and dense (right lower corner) and become hot and during the black hole flyby. If 
the time scale on which the white dwarf can react (dynamical time scale Iq) is longer than 
the burning time scale (tb, see Eq. f|T8l) ). the star cannot expand rapidly enough to quench 
burning. 
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Fig. 15. — Energy generated in nuclear burning (in units of the star's binding energy) as 
a function of time (same simulations as referred to in Figures [13] and [T^ . Here time is 
measured in units of the dynamical timescale (to) of the initial 0.2 Mq white dwarf. In 
addition, the time axis has been shifted so that the maximum E\, occurs at t = 0. 
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Fig. 16. — A 1.2 Mq white dwarf (modeled with more than 5 x 10^ SPH particles) approaching 
a 500 Mq black hole on a parabolic orbit with pericenter distance rmin = rT/3.2. Shown are 
temperature (in units of 10^ K) and surface density (in cgs units) cuts at various instants 
along the orbital (xy-) plane. 
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Fig. 17. — Evolution of the abundances during the disruption and ignition of a 1.2 Mq white 
dwarf approaching a SOOMq black hole (same simulation as referred to in Figure [T6|) . 
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Fig. 18. — Comparison of the average properties of the hottest 10% of the SPH particles for 
an initially cold and warm star (runs 8 and 9 in Table 1). 
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Fig. 19. — The evolution of the compressed, and tidally disrupted white dwarfs of various 
masses in the p — T plane: A 0.2 He Mq white dwarf approaching a lO^M© black hole with 
/3 = 12; a 0.6 C/0 Mq white dwarf approaching a 10'^ Mq black hole with j3 = 1.5 and a 
1.2 C/0 Mq white dwarf approaching a lO^M© black hole with /? = 1.5. The hottest 10% 
of the particles are identified and their average temperature (in units of 10^ K) is plotted 
as a function of their average density (in cgs units). These trajectories always start cold 
and dense (right lower corner) and become hot and during the black hole flyby. For the 1.2 
Mq and 0.6 Mq cases, the time scale on which the white dwarf can react is similar to the 
burning time scale and the star expands rapidly enough to quench explosive energy release. 
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Fig. 20. — Frequency and amplitude of gravitational radiation produced by the disruption 
of white dwarfs approaching an intermediate mass black hole, assuming (3 — 1 and that the 
cluster hosting the black hole is at D = 8kpc. 
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Fig. 21. — Gravitational wave amplitudes. Top panel: gravitational radiation produced 
by the disruption of a 0.2 Mq white dwarf approaching a lO^M© black hole with /3 = 12. 
Bottom panel: gravitational radiation produced by the disruption of a 0.6 Mq white dwarf 
approaching a lO^M© black hole with (3 = 1.5. 



